Alternative quantization of the Hamiltonian in isotropic loop quantum cosmology 



ON 

o 
o 

(N 

< 

oo 
(N 



o 

i 

wo. 



> 

m 

On 

(N 
O 
On 
O 



X 



Jinsong Yang 1 , You Ding 1 ' 2 , and Yongge MaQ 
1 Department of Physics, Beijing Normal University, Beijing 100875, China 
2 Centre de Physique Theorique de Luminy, Universite de la Mediterranee, F-13288 Marseille, EU 

Since there are quantization ambiguities in constructing the Hamiltonian constraint operator in 
isotropic loop quantum cosmology, it is crucial to check whether the key features of loop quantum 
cosmology, such as the quantum bounce and effective scenario, are robust against the ambiguities. 
In this paper, we consider a typical quantization ambiguity arising from the quantization of the 
field strength of the gravitational connection. An alternative Hamiltonian constraint operator is 
constructed, which is shown to have the correct classical limit by the semiclassical analysis. The ef- 
fective Hamiltonian incorporating higher order quantum corrections is also obtained. In the spatially 
flat FRW model with a massless scalar field, the classical big bang is again replaced by a quantum 
bounce. Moreover, there are still great possibilities for the expanding universe to recollapse due to 
the quantum gravity effect. Thus, these key features are robust against this quantization ambiguity. 
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I. INTRODUCTION 

An important motivation of the theoretical search for a 
quantum theory of gravity is the expectation that the sin- 
gularities predicted by classical general relativity would 
be resolved by the quantum gravity theory. This expec- 
tation has been confirmed by the recent study of cer- 
tain isotropic models in loop quantum cosmology (LQC) 
[3, [2, El, which is a simplified symmetry-reduced model of 
a full background-independent quantum theoryofgravity 
[U , known as loop quantum gravity (LQG) . In 

loop quantum cosmological scenario for a universe filled 
with a massless scalar field, the classical singularity gets 
replaced by a quantum bounce [!, [H, [nj. Moreover, it 
is revealed in the effective scenarios that there are great 
possibilities for a spatially fiat FRW expanding universe 
to recollapse due to the quantum gravity effect QJJ • How- 
ever, as in the ordinary quantization procedure, there 
are quantization ambiguities in constructing the Hamilto- 
nian constraint operator. Thus a crucial question arises. 
Whether the above significant results come from certain 
particular treatment of quantization? Before confirming 
the robustness of the key results against the quantization 
ambiguities, one could not believe that they are not some 
artifact under particular assumptions. 

In this paper, we consider a typical quantization ambi- 
guity arising from the quantization of the field strength 
of the gravitational connection. An alternative Hamil- 
tonian constraint operator is constructed. Semiclassical 
states are then employed to show that the new Hamil- 
tonian operator has the correct classical limit. The ef- 
fective Hamiltonian incorporating higher order quantum 
corrections is also obtained. In the spatially flat FRW 
model with a massless scalar field, the classical big bang 
is again replaced by a quantum bounce. Moreover, there 
are still great possibilities for the expanding universe to 



recollapse due to the quantum gravity effect. 

In the spatially fiat isotropic model of LQC, one has 
to first introduce an elementary cell V and restrict all 
integrations to this cell. Fix a fiducial flat metric °q ab 
and denote by V the volume of the elementary cell V in 
this geometry. The gravitational phase space variables — 
the connections A l a and the density- weighted triads Ef — 
can be expressed as 



At 



cV-^ 3o u l a and E?=pV~ 2 / 3 



°q °e?, 



(1) 



where (°ui' l a , ef) are a set of orthonormal co-triads and 
triads compatible with °q ab and adapted to the edges of 
the elementary cell V. The basic (nonvanishing) Poisson 
bracket is given by 



{c, p} 



8ttG7 



(2) 



where G is the Newton's constant and 7 is the Barbero- 
Immirzi parameter. 

To pass to the quantum theory, one constructs a kine- 
matical Hilbert space Tif^ = L 2 (MBohr, d^Bohr), where 
KBohr is the Bohr compactification of the real line and 
d/^Bohr is the Haar measure on it The abstract *- 

algebra represented on the Hilbert space is based on the 
holonomies of connection A l a . In the Hamiltonian con- 
straint of LQG, the gravitational connection A l a appears 
through its curvature F* b . Since there exists no operator 
corresponding to c, only holonomy operators are well de- 
fined. Hence one is led to express the curvature in terms 
of holonomies. Similarly, in the improved dynamics set- 
ting of LQC [3j, to express the curvature one employed 
the holonomies 



1) := cos — 11 + 2 sin — n 
2 2 



(3) 
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along an edge parallel to the triad °ef of length fiW \p\ 
with respect to the physical metric q ab , where I is the 
identity 2x2 matrix and ty = — i<Ji/2 (o"j are the Pauli 



2 



matrices). Thus, the elementary variables could be taken 
as the functions exp(i/Ic/2) and the physical volume 
V = |p| 3 / 2 of the cell, which have unambiguous opera- 
tor analogs. 



II. THE ALTERNATIVE HAMILTONIAN 
CONSTRAINT OPERATOR 

In general relativity, the dynamics of a gravitational 
system is determined by the Hamiltonian constraint. 
Many problems, such as the big-bang singulary, arise in 
classical dynamics. One expects that some quantum dy- 
namics can resolve these problems. Hence, it is important 
to have a well-defined Hamiltonian constraint operator in 
LQC. An improved Hamiltonian constraint operator has 
been constructed in However, there are quantization 
ambiguities in the construction. In this section, we will 
construct an alternative Hamiltonian constraint operator 
by a different quantization procedure. 

Because of spatial flatness and homogeneity, the gravi- 
tational part of the Hamiltonian constraint of full general 
relativity is simplified to the form 



Cgrav — 
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d 3 xNe ijk F* b e- 1 E^E bk , 



(4) 



where e := y/\ det_E|, the lapse N is constant and we will 
set it to one. 

The procedure used in LQC (and specifically in [3LIT2I]) 
so far can be summarized as follows. The term involving 
the triads can be written as 



E aj E' 



bk 



sgn(p) 



x Tr (hf {h 



1/3 



^abc o ,k 
6 UJ^ 



V}ri). (5) 



To express the curvature components F^ b in terms of 
holonomies, one considers a square in the i-j plane 
spanned by a face of V, each of whose sides has length 
M\/[pI with respect to q a b- Then the ab component of the 
curvature is given by 



ib'i 



lim 

Am- 



1 



>0 ft 2 Vo 



2/3 



°, i l °, J 



(6) 



where Ar a is the area of the square under consideration, 
and the holonomy around the square □ y is just the 
product of holonomies along the four edges of □ y , 



(7) 



However, quantization ambiguities arise here, since the 
approach to express the curvature components F* b in 
terms of holonomies is not unique. Hence the corre- 
sponding operators in different approaches will be dif- 
ferent from each other. In the following, we will con- 



sider an expression of the curvature components differ- 
ent from Eq. ^ . Taking account of the definition of 
holonomies, we have the identity 



(fi)- 1 



lim ■ 

/2->0 



/' 



= lim 4sin(/2G/2)T ' = 2cn. (8) 



Hence the curvature of connection can be written in 
terms of the holomomies as 



lim 



» 



2fi 2 Vo' 



U [a ^b]- 



(9) 



Combining Eqs. {5]) and (J9j> , the Hamiltonian constraint 
can be written as 



CWav — lim 



sgn(p) 



. e ijk Tl 



- h 



lim C R 



(10) 



Since the constraint is now expressed in terms of ele- 
mentary variables and their Poisson bracket, it can be 
promoted to a quantum operator directly. The resulting 
alternative regulated constraint operator with symmetric 
factor-ordering reads 



c r 



(P) 



fIC 



12isgn(p) / . fxc y fic 
'sin— V cos — 



i 3 ei v 



cos — V sin — 
2 2 



sin ■ 



(11) 



where, for clarity, we have suppressed hats over the oper- 



and i\ = VGh. 



To 



ators h\ , sin(/xc/2) and cos(/2c/2) 

deal with the regulator ft, we adopt the improved scheme 
0| . We shrink the length of holonomy edges, as mea- 
sured by the physical metric q a b, to the value \/A, where 
A = 4v / 37T7^p is a minimum nonzero eigenvalue of the 
area operator [lb] ]. Thus we are led to choose for p, a 
specific function p,(p), given by 



(12) 



It is convenient to work with the w-representation. In 
this representation, states \v) constituting an orthonor- 
mal basis in 7ij^ v is more directly adapted to the volume 
operator V, 



V\v) 



Stt 7 £ 



3/2 



K 



\v), 



(13) 
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where 



K 



3A ' 



The action of exp(z/j,c/2) is given by 

exp(ip,c/2)\v) = \v + 1) 



(14) 



(15) 



Hence the alternative Hamiltonian constraint operator is 
given by 



r R -sin^ 



12i sgn(w) 



sin — K cos — 
V 2 2 



cos — V sin — 
2 2 



(AC 



=: sin — AA sm — . 
2 2 



(16) 



It is easy to show that A is well defined and \v) is an 
eigenvector of A. Furthermore, the eigenvalues of A are 
real and negative. So A is a negative definite self-adjoint 
operator on Wj^ v . Hence, C^ lav is a negative-definite 

self-adjoint operator on . The action of C^. m on 

the basis |u) of Wf[^ v is given by 

Cg R ra » = /»|« + 2) + f {v)\v) + f_{v)\v - 2), (17) 
where 

27 l8irK£ v 



\(v + l)(\v 



H). 



+ 4 V 6 7 3 / 2 ' 

/i(«) = /> -2), /»=-/» -f'_(v). (18) 

Thus, C R av is again a difference operator. Recall that by 
contrast to Eq. (|16|) . the Hamiltonian constraint opera- 
tor defined in [3j reads 



sin(/ic) 



3z sgn(u) 



7T7 3 /X 3 £ 2 



cos — V sm — 
2 2 



sm — V cos — 
V 2 2 

sin(/2c) 



sin(/ic)A sin(/ic). 



(19) 



This shows a quantization ambiguity arising from the 
quantization of the field strength of the gravitational con- 
nection. 



To identify a dynamical matter field as an internal 
clock, we take a massless scalar field <j> with Hamiltonian 
Cgj, — |p|~ 3 / 2 p|/2, where denotes the momentum of 
0. In the standard Schrodinger representation, the mat- 
ter part of the quantum Hamiltonian constraint reads 

CV> = |p| — 3 / 2 /2. Thus we get the total constraint as 



C R 



1 

leTrG^grav 



III. THE CLASSICAL LIMIT AND THE 
MODIFIED FRIEDMANN EQUATION 



It has been shown in [TJ [l3| that the improved Hamil- 
tonian constraint operator constructed in [3j has the cor- 
rect classical limit. In this section, we will show that the 
alternative Hamiltonian constraint operator constructed 
in last section also has the correct classical limit. More- 
over, the effective Hamiltonian incorporating higher or- 
der quantum corrections can also be obtained. In order to 
do the semiclassical analysis, it is convenient to introduce 
new conjugate variables by a canonical transformation of 
(c,p) as 



and 



sgn(p)|p| 3/2 



(20) 



with the Poisson bracket {b, v} = 1/H. In terms of these 
new variables, the classical Hamiltonian constraint can 
be written as 



C 




Let us first consider the gravitational part. Since there 
are uncountable basis vectors, the natural Gaussian semi- 
classical states live in the algebraic dual space of some 
dense set in Hj^ v . A semiclassical state {^(b ,v ) I peaked 
at a point (b ,v ) of the gravitational classical phase 
space reads: 



(*(6»,t 



-[(«-• v ) 2 /2d 2 ] ibo(v-v ) 



01= (22) 



where d is the characteristic "width" of the coherent 
state. For practical calculations, we use the shadow of 
the semiclassical state (^(b ,v a )\ on the regular lattice 
with spacing 1 [14j , which is given by 



i*) = E [■ 

riGZ 



-{e 2 /2)(n-N) 2 -i(n-N)b 



|n + A), (23) 



where A € [0, 1), e = 1/d and we choose v = N + A, 
here N G Z. Since we consider large volumes and late 
times, the relative quantum fluctuations in the volume of 
the universe must be very small. Therefore we have the 
restrictions: 1/N <C e <C 1 and 6„ < 1. One can check 
that the state (f22|) is sharply peaked at (b „, v„ ) and the 
fluctuations are within specified tolerance [ll|,ll3|. The 
semiclassical state of matter part is given by the standard 
coherent state 



)e -[(</>-<M 2 /2<r 2 ] (24) 
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where a is the width of the Gaussian. Thus the whole 
semiclassical state reads {^(b ,v )\ ®(^ r W <) ,zv)l- 

The task is to use this semiclassical state to calcu- 
late the expectation value of the Hamiltonian operator 
to a certain accuracy. In the calculation of (C^ av ), one 
gets the expression with the absolute values, which is 
not analytical. To overcome the difficulty we separate 
the expression into a sum of two terms: one is analyti- 
cal and hence can be calculated straightforwardly, while 
the other becomes exponentially decayed out. We thus 
obtain (see the Appendix for details) 



(Cgrav) — — 54 



6 7 3/2 ^°' 



e e2 sin 2 b + - ( I - < 



(25) 



In the calculation of (C$), one has to calculate the expec- 
tation value of the operator |p|~ 3 / 2 . A straightforward 
calculation gives: 



-3/2) 



/' 6 



3/2 



K r 



\v \i 2MV 9M 2 

+ 0(1/H 4 e 4 )] + 0(e- N °<°) +0(e-^). 

(26) 

Collecting these results we can express the expectation 
value of the total Hamiltonian constraint, up to correc- 
tions of order l/|i> | 4 e 4 and e" 7 ' / e , as follows: 



(c R )=- 




sin b Q 



3/2 



K 



-(—) ~ 

2 \ton&) K\ 



2a 2 



1 



2\v n \ 2 e 2 



(27) 



Hence the classical constraint (|2"Tj) is reproduced up to 
small quantum corrections. Therefore, the new Hamil- 
tonian operator is also a viable quantization of the clas- 
sical expression. For clarity, we will suppress the label 
o in the following. Using the expectation value of the 
Hamiltonian operator in Eq. (|27p . we can further ob- 
tain an effective Hamiltonian with the relevant quantum 
corrections of order e 2 , l/v 2 e 2 7 h 2 /a 2 as 



27 8k Kl 



n ° s ' 8kG V 6 7 3 / 2 



v\ I sin 2 b H — e 2 



'87r 7 ^ 2 \ 3/2 | v | h 2 



K P y 1 ' 2a 2 p 2 ' 2v 2 e 2 J ' 

(28) 



where p—\ ( 8 J^ e i ) (if)' * P% is the density of the mat- 



ter field. Then we obtain the Hamiltonian evolution 
equation for v by taking its Poisson bracket with H^ s 
as 



v = {v,H%} = 



27K 



2^Gtvf 
Further, the vanishing of Tl^g implies 

h 2 1 



v \ sin b cos b. 



(29) 



sin* 6=^-1 



(30) 



p' c V 2cr 2 p 2 2v 2 e 2 J 2 
where p' c — 3 / (2ttGj 2 A) . The modified Friedmann equa- 



tion can then be derived from Eqs. 

2 



and ([3H)l as 




(31) 



Recall that, up to the quantum fluctuation of matter 
field, the modified Friedmann equation given in (llj reads 



H 2 = 



InG 



1-^11 

Pc 



1 



1 



2v 2 e 2 



-2e 2 ^ 



(32) 



where p c = 3/(8ttG7 2 A). Comparing Eq. ((3T]) with 
([32|) . we find that the leading order critical energy density 
reads p' cl . it — 4p cr i t . We can also express the new modified 
Friedmann equation by using p c as 



H 2 




(33) 



IV. DISCUSSION 

Because the properties of the alternative Hamiltonian 
constraint operator in Eq. (|16p are similar to the one 
in Eq. (Tj"9|) . the physical Hilbert space, Dirac observ- 
ables and so on investigated in Q can also be straight- 
forwardly obtained. Although there are quantitative dif- 
ferences between the two versions of quantum dynam- 
ics, qualitatively they have the same dynamical features. 
In the leading order approximation, the universe would 
bounce again from the contracting branch to the expand- 
ing branch when the energy density of scalar field reaches 
to the critical p' c ^ t = 4p cr ;t. The quantum bounce implied 
by Eq. (|3Tj) is shown in Fig. [1] 

On the other hand, the key results discussed in [ll| 
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1,000 2,000 3,000 4,000 5.000 6.000 7,000 8,000 9,000 10,000 



FIG. 1: The effective dynamics represented by the observable 
v\ t j > are compared to classical trajectories. In this simulation, 
the parameters were: G — h — 1 , 
0.01 with initial data v a = 100 000 . 



10 000, e = 0.001, a = 



can be carried out similarly. It is easy to see from Eq. 
(I3ip that, while the term containing the quantum fluc- 
tuations of matter field is qualitatively negligible, the 
asymptotic behavior of the quantum geometric fluctu- 
ations plays a key role for the fate of the universe. By 
the ansatz e = a(r)v~ r (^ with < r{(j>) < 1, there are 
great possibilities for the expanding universe to undergo 
a recollape in the future. The recollape can happen pro- 
vided < r < 1 in the large scale limit. Suppose that the 
semiclassicality of our coherent state is maintained asym- 
totically. This means that the quantum fluctuation l/e 
of v cannot increase as v unboundedly as v approaches 
infinity. Thus the recollape is in all probability as viewed 
from the parameter space of r (</>). For example, in the 
scenario when r = asymtotically, besides the quan- 
tum bounce when the matter density p increases to the 
Planck scale, the universe would also undergo a recol- 
lapse when p decreases to p co ii ~ e 2 p' c /2. Therefore, the 
quantum fluctuations again lead to a cyclic universe in 
this case. The cyclic universe in this effective scenario 
is illustrated in Fig. O This is an amazing possibility 
that quantum gravity manifests herself in the large scale 
cosmology. Nevertheless, the condition that the semi- 
classicality is maintained in the large scale limit has not 
been confirmed. Hence further numerical and analytic 
investigations to the properties of dynamical semiclassi- 
cal states in the model are still desirable. It should be 
noted that in some simplified completely solvable models 
of LQC (see @ and [IH), the dynamical coherent states 
could be obtained, where r(<f>) approaches 1 in the large 
scale limit. While those treatments lead to the quantum 
dynamics different from ours, they raise caveats to the 
inferred re-collapse. 



FIG. 2: The cyclic model is compared with expanding and 
contracting classical trajectories. In this simulation, the pa- 
rameters were: G = h = 1 , p = 10 000 , e = 0.001, a = 0.01 
with initial data v — 100 000. 



In conclusion, the key features of LQC in this model, 
that the big bang singularity is replaced by a quantum 
bounce and there are great possibilities for an expanding 
universe to recollapse, are robust against the quantiza- 
tion ambiguity which we have considered. 



ACKNOWLEDGMENTS 

This work is a part of project 10675019 supported by 
NSFC. 



APPENDIX 

Let us calculate the expectation value of C^. 



grav i 



<3 



R \ 
rav / 



<*!*> 



(34) 



Applying the Poisson resummation formula 

OO OO n QQ 

]T g(n + x)= e2 ™ X / 9(y)e- 2 ™ y dy (35) 

11 — — OO 71 — — OO — OO 

to the norm of the shadow state (f!?3")) . one obtains 

(*|*> =^e- e2 " 2 = ^ (l + 0(e - £)) . (36) 

By Eq. (|17[) . we obtain the action of the gravitational 
Hamiltonian operator on the shadow state as 



C' g R r avl*)=E 



-(e 2 /2)(n-N) 2 -i(n-N)b 
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x f + (n + A)|n + A + 2) + f' {n + A)|n + A) 
+ /i(n + A)|n + A-2) . 
Then a straightforward calculation shows that 

(*ia R rav l*> 



\n-NY 
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4 V 6 7 3 / 2 



2e" e cos(26 )/|(n + A - 1) 

-(/;(n + A) + /i(n + A))" 
2e" e2 cos(26 o )5 - (5i + , (37) 



where 



/ > ) = fvf?7l (n + 1)(|n + 2| - H) ' 

f'-^)- 2 i^^(n~l){\n\-\n-2\), 
and we have set 

x (|n + A + m + 1| - \n + A + m - 1|), 
S m :=2 ^ e^™^) 2 ( n + A + to) 

=2(iV + A + to)^^' 2 ™ 2 - 



We may bound |<!5 m ,fc — S m \ by an exponentially sup- 
pressed term: 



I S m S m | — 



net 



< 



^ e - e 2 („-A0 2 (n + A + m) 

x (|n + A + to + 1| 
- |n + A + m- 1| - 2) 



5] e-^-^V + A + m) 

-Kn+A+m<l 

x (2(n + A + to) - 2) 

+ 4 51 e _f2(n_Ar)2 (n + A + to) 

n+A+m< — 1 



53 e- e2 ("-^ 2 (n + A + TO) 

n+m,— — 1 



x (2(n + A + m) - 2) 



+ 4 



E 



n+ m< — 1 



^ E 



,-e 2 (n-JV)' 



\n-NY 



(n + m) 



n+m<— 1 

e- e2 ( m+JV ) 2 A(A-l) 

+ e- e2 ( m+JV ) 2 (A-l)(A-2) 



453 P • J i»- -» - x ) 
?i>i 



+ 4A53 



e v ' n 

t 2 (n+m+N) 2 



(38) 



Using the Euler-Maclaurin summation, for some positive 
constant a independent of N, m or e, we obtain: 

/>oo 

4 53 e - e2 ("+ m + w ) 2 n ^ae-^ 2 + 4 / e^+^+^xdx 



<ae 



-N- 



e 3 (m + AT) 



2 (m+iV) 2 



(39) 



Hence we obtain \S m — S m \ = 0(e e2jv2 ). Therefore, we 
have 



2e _e cos(26 o )S 



4 V 6 7 3 / 2 
-(Si + S-ijj+OCe-^). (40) 

Finally, we can collect the above terms to calculate the 
expectation value (J34j). Using (*|*) = E„ eZ e - ^" 5 \ we 
have 



(C g R rav ) = 27^^ (N + A) [e- a cos(26 G ) - 1 

+ o( e - w2e2 ) 



6 ^ 11/01 
0( e ^ 2 ). 



e" e sin 2 6 + -(l-e- £ ) 
(41) 
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